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Abstract 

We study the critical behaviour near the threshold where a first bound state 
appears at some value of coupling constant in an attractive short-range poten- 
tial in 2 + e dimensions. We obtain general expression for the binding energy 
near the threshold and also demonstrate that the critical region is correctly 
described by an effective separable potential. The critical exponent of the 
radius of weakly bound state is shown to coincide with the correlation length 
exponent for the spin model in the large- iV limit. In two dimensions, where 
the binding energy is exponentially small in coupling constant, we obtain a 
general analytic expression for the prefactor. 
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I. INTRODUCTION 



It is well known, that in three dimensions a bound state for a particle in a short-range 
potential exists not at an arbitrary value of the coupling constant A, but only at A > A c , 
where A c is a critical value which depends on the particular potential. This may be viewed 
as the simplest example of a quantum phase transition, when e.g. an excitation gap vanishes 
as some parameter of the Hamiltonian is varied (see e.g. Q for a review). In a sense, such a 
behaviour is similar to the second-order phase transition. Near the threshold the energy of 
the bound state behaves like an 'order parameter' E ~ (A — A c )^ , where A plays the role of 
temperature and (3 is the critical exponent. Deeper investigation of this critical behaviour 
is interesting in itself and may have some applications (see || and references therein). 

The lower critical dimensionality for this transition is d — 2, since in two dimensions 
there always exists a bound state with energy exponentially small in A ||, hence A c = in 
this case. For this reason it seems natural to study the critical behaviour in 2 + e dimensions, 
considering e to be a small parameter, as it was done for the phase transition in the nonlinear 
O(N) sigma model || and also in the theory of Anderson localization (see e.g. ||). It appears 
possible to develop an e-expansion both for the wavefunction at the critical point and for 
the critical coupling A c (e) Even the first two terms of the expansion of A c (e) in powers 

of e give rather accurate estimate for A c in three dimensions. 

In this note we consider the onset of the first bound state in 2 + e dimensions in more 
detail. First, we demonstrate that (3 = 2/e for e < 2 and f3 = 1 above four dimensions. The 
result at small e is not unexpected, since (3 must go to infinity as e — ► to reproduce the 
exponential dependence of E on A in two dimensions. This expressions for j3 are consistent 
with the results of Lassaut et al [Q]. Though they have studied three dimensional case with 
non-zero orbital momentum I, this is equivalent to the s-state problem in 2 + e dimensions 
with e = 1 + 21 (see below) . 

However, the method used here is different from that of 0. Starting from the integral 
representation of the Schrodinger equation we first derive exact results for the binding energy 
at A ~ A c . Then we show also that the correct description of the critical region is given by the 
separable approximation, with the true interaction potential replaced by an effective nonlocal 
separable one, which depends on the zero energy solution of the Schrodinger equation at 
A = A C . 

When A —>■ A c the radius of the bound state diverges as (A — \ c )~ u with the critical 
exponent v — Therefore the result obtained may be represented in terms of z/, namely 
v = 1/eate < 2 and v = | at e > 2. This exponent coincides with the correlation 
length exponent in the spherical model (equivalent to the iV-component spin model at iV — > 
oo) || and with the localization length exponent in the self-consistent theory of Anderson 
localization |J. It is not clear whether this coincidence imply some nontrivial relation 
between these models, but still it seems to be rather interesting. 

The same critical exponents were obtained also from scaling considerations by Hwa and 



Nattermann [[KJ and Kolomeisky and Straley [IT], who considered the problem of unbinding 
a directed polymer from a columnar defect in the presence of quenched disorder. Since such 
a polymer may be viewed as a worldline of a quantum particle, in the clean case this problem 
is essentially the same as the one discussed here. 

Next, the approach used here also makes it possible to obtain an asymptotic expression 
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for the energy of the bound state in two dimensions at A — > along with the preexponential 
factor. This general analytic expression for the prefactor seems, to the best of our knowledge, 
to be a new one. We also calculate this prefactor for some simple potentials and discuss its 
connection with the e-expansion for A c . 

As a direct application of the general result for the binding energy in two dimensions 
we also consider a case of a two-center potential. The attractive force between the centers 
due to the bound state is shown to be of the Coulomb type (compare with ]T2||), and the 
universal prefactor is found without solving the Schrodinger equation. 



II. CRITICAL BEHAVIOUR IN 2 + e DIMENSIONS. 

Consider the Schrodinger equation in d dimensions 

- Atf(r) +AV(r)#(r) = £tf(r), (1) 

( h — 1, 2m = 1), where A is the coupling constant and V(r) is a short-range attractive 
potential. We assume V(r) to decrease faster then 1/r 2 as r — > oo and shall deal here with 
the weakly bound state of the size ~ ( y —E)~ 1 / 2 ^> a, where a is the radius of the potential. 

The first bound state in the problem appears at some critical value of the coupling 
constant A = A c . Let us denote by ipo the wavefunction of this state at the threshold. Then 
ipo obeys the following zero energy equation 

- A^o + A c ^(r)^ = (2) 
It is convenient to normalize i/jq by the condition 

j dv^{v)V{v) = -l (3) 

Note, that ipo(r) need not be square integrable and the convergence of the normalization 
integral in (|3]) is garantied by the short-range potential V(r). 

We shall now determine both the critical value A c and the binding energy at A ~ A c . For 
this purpose first rewrite the original Schrodinger equation in the integral form 

tf( r ) = -A J dr'G E {r - r')V(r')V(r') , (4) 

where Ge{t) is the Green's function of the free particle. This integral equation may be 
viewed also as an equation determining X(E) for the given negative energy E of the bound 
state. Multiplying both sides of by V(r)\l/(r) and integrating over r we easily obtain 

I J dr dr' V(t)V(t) G E (r - r') V (r')^(r') 

A~ Jdr^ 2 (r)V(r) ^ ' 

If we put here E = then ^ — > ipo and using the normalizing condition (|3]), we have 

— = [ drdr'W(r)G (r-r')W(r'), (6) 
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where 



W(r) = V(r)Vo(r) . (7) 
At zero energy Gq(t) is merely the Green function of the Laplace operator and 

CM = ~ (8) 

(see e.g. fll3f), where 

27T 1+e / 2 



Or 



r(i + e /2) 



(9) 



is the area of the unit sphere in 2 + e dimensions. Then the critical value of the coupling 
constant may be represented as follows 



A, 



27T 1+e / 2 

r(l + e/2) 



dr dr 



f W(r) W(r') 



r — r 



(10) 



This expresion explicitly demonstrate, that normally A c tends to zero as e when we approach 
two dimensions. 

Now, in the vicinity of the critical point we may write 



# = -00 + 5V , G E = G + 5G E 



(11) 



and assume all corrections to zero energy values to be small. If we substitute ( [LTD i n 
and retain only terms of first order in 5^ we finally obtain 



_ = _+/" drdr'W(r)5G E (r - r')W(r') 

A \ r J 



(12) 



with 1/A C given by @. Note that the terms containing 8^ cancel out. This cancellation is 
a consequence of the zero energy equation 



^o(r) = -X C J dr'G (v - r')V{r')M*') 



(13) 



and is actually due to the right-hand side of ([5]) being a variational functional which is stable 



against small variations of the true wavefunction [IJJ]. In the limit E — > we can expand 
5Ge(y — r') in equation ( |12D in powers of E. In 2 + e dimensions with e < 2 we have at 
small negative E 



6G, 



rq-e/2) (Ey> 

2evr \ Aix) V ' 



(14) 



(see Appendix A). Substituting this expression in ( J12|) , we finally obtain 



E = — I A 



X-Xr 



2/6 



< e < 2 



(15) 



where 



Til + e/2) jdvdv'W{v)W{v')\v-v'\- e nfi , 
r(l-e/2) (JdrW(r)) 2 ' 1 j 

At e > 2 , the leading term in the expansion fll4]) for 5Ge{t) is ~ E , and 

£~(A-A C ) e>2. (17) 

These equations formally solve the problem of the critical behavior near the transition where 
the first bound state appears in d = 2 + e dimensions. In three dimensions equation ([15]) 
leads to E ~ (A — A c ) 2 and for the square well potential one can easily verify (using ip from 
equation (P^)of Appendix B) that equations ([TDD , ( P^D gi ye the correct answer A c = 7r 2 /4a 2 
and A = it 2 /8a. Note, that if the Schrodinger equation is solved in the critical point, i.e. A c 
and ipQ are known, one can also evaluate the prefactor A. 

The results obtained should be compared with that of 0, where the /-wave case was 
considered, because the radial s-wave Schrodinger equation in d = 2 + e dimensions is 
equivalent to the three dimensional equation with non-zero orbital moment / = (e — l)/2 
(see Appendix A). Hence e.g. the dependence (P~5[) is the same, as E ~ (A — \ c ) 2 ^ 2l+l ^ 
obtained in pf for I < 1/2. 



The particular form of equation ( I4) suggests that the correct description of the critical 
region near the threshold A ~ A c can be obtained within the separable approximation. This 
approximation, widely used in nuclear physics, involves replacing the original potential V(r) 
with a non-local separable one, for which the Schrodinger equation is exactly solvable. In 
our case one should take 

V sep = -V\if> )(i> \V, (18) 

where ipo is the zero energy solution normalized by the condition (|3|). If the ground state 
wave function \I/ « ip then V sep is in a sense close to V, since (V — V sep )il)o = 0. 
The Schrodinger equation for a particle in the potential QI8"D reads 



- Atf(r) - AV(r)^ (r) J drV(r>„(r')f (r') = EV(r) (19) 

and has an obvious solution for the bound state, which up to a normalizing constant is given 
by 

#(r) = -A J dr'G E {r-v')W{v') , (20) 

The energy of the bound state is determined by substitution of ( |2"0"D into equation ( |T9"D i.e. 
from the equation 

1 = A J dr dr' W[r) G E (r - r') W{r r ) . (21) 

This is just the same equation as (|T2"|), since G E = G + 5G E and 1/A C is determined from 
(|). Therefore the separable approximation ([TJ) results in exact expressions (|15]) , (filf) in 
the close vicinity of the critical point. The validity of the separable approximation seems 



5 



to be due to the wavefunction ( p0|) having correct asymptotic behaviour at r > a. This 
is similar to the one dimensional case, where the energy of the weakly bound state can be 
obtained replacing the true V(x) with a suitable 5-function potential, which also may be 
viewed as a separable one. In Appendix B we show how one can naturally arrive at the 
separable potential of the form ([Op. 

As A — ► A c the radius of the bound state £ ~ (— E)^ 1 ! 2 goes to infinity as (A — \ c )~ u , 
where we have introduced a new critical exponent v. Then from ( |T5D and flTT| ) it follows 
that 



T-2 2<^<4 



(22) 



Critical exponent v diverges as d — > 2 and 'freezes' above d = 4 at the mean field value 
= ~. This is precisely the correlation length exponent for the iV-component spin model 
at N — > oo 0. Another model where equation (^2j) arises, is the self-consistent theory 
of Anderson localization for a particle in a random potential. In this case the localization 
length diverges as (Ep — E c )~ v if the Fermi energy Ep approaches the mobility edge E c and 
v is also given by equation (|22"D . It was even argued that this result for v is valid beyond 
the self-consistent approximation and is in fact an exact one |l5|Jl6 l. 



This interesting coincidence arises from the fact that in all these models resulting equa- 
tion, determining the behaviour of the correlation length £, has the form similar to (|21~|) with 
E — > C,~ 2 . For example, in the self-consistent theory of Anderson localization the localization 
length £ is given by the equation 

l = BE ?~r«gL, (23) 

Jo q + 

where B is some constant and go is a momentum cut-off |J. Comparing ( f23| ) with ([H|) we 
see that this is indeed the equation for the binding energy in an effective short-range 
separable potential with A ~ 1/Ef^ e , a ~ l/q Q . If the Fermi energy increases, A tends to 
zero and at A = A c the bound state disappears. This critical point obviously corresponds to 
the Anderson transition. Perhaps this is not a mere coincidence and some direct mapping 
between these models might be established. 

To conclude this section we should like to mention that at small e the critical exponents 
derived here can be obtained without actually solving the Schrodinger equation. This was 



done e.g. by Hwa and Nattermann JlIJ and Kolomeisky and Straley [[Ll|| , who considered 
the problem of unbinding a directed polimer (i.e. the worldline of a quantum particle) from 
a columnar defect. In this case simple scaling arguments immediately lead to equation (J22J). 
In fact, one can take any quantity (not necessarily the free energy as in pO|JTT|| ), depending 
on some scale and look at the perturbation theory in A. Consider e.g. the Born series for 
the s-wave scattering amplitude f(k) (see e.g. JTj 



f(k) = f 1 (k)+f 2 (k) + --- (24) 

where at small k 

/i ~ XV(0) , W'jf^J^L. (25) 
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and V(q) is the Fourier transform of the potential. For the perturbation theory to be valid 
it is necessary that /2//1 «C 1. At small e one has from (p5|) / 2 ~ A 2 t^ 2 (0) (1 — (k/k ) e )/e, 
where ko ~ 1/a, and hence the particle is essentially free on a scale A; if 



k 
h 



A 



Uo, 



< 1 



(26) 



(compare with |10|), where A c ~ eV(0). If A > A c then (|26|) is not fulfilled at k = and the 
particle is bound. But on a length scale 



k^i 



-1 



AY) 



'A-Ac 



(27) 



one can neglect the potential. Hence £ from ( 27D may be viewed as the radius of the bound 
state. The critical exponent obtained in this way is the same as (|22|). 



III. WEAKLY BOUND STATES IN TWO DIMENSIONS. 

In three dimensions neither A c nor A are known exactly for an arbitrary potential, since 
we can solve the zero-energy problem only in some special cases. But the situation is different 
in two dimensions. In this case arbitrarily weak short-range attractive potential binds a 
particle, so that A c = 0. Then the solution to zero-energy equation (fj) is obviously 4>o(r) = 
const, which, according to the normalization condition (|j), results in ij)Q = I / drV^)]^ 1 ^ 2 
and 

w ^ = \i*vl)\* (28) 

Next, in the limit e->0we have A c ~ e and A — > 1. Therefore in this limit the right- 



hand side of equation (|lg) for the binding energy turns to an exponential function. Then 
for k = \J —E we have 

^ c ^(-w^m)' (29) 

where the prefactor C is determined from the expansion of A in powers of e 

A~ l + elnC + --- . (30) 
Expanding A from fllB]) in e we obtain 



/ J drdr'V (r)V(r') ln|r-r'|\ 

C = exp r 7 U*v<T)y -) ' (31) 

where 7 = 0.577. . . is the Euler's constant. Thus, in two dimensions we have a general 
explicit expression for the energy of the weakly bound state. In contrast to the one dimen- 
sional case, where at a —>■ we may approximate any short range potential by the 5-function 
and the binding energy depends only on one potential dependent integral / dxV(x), here 
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TABLE I. Values for Ai and A2 from (|32]), a from (|36| ) and the prefactor C from (|3l| ) for 
several different potentials V(r). 



-y(r) 


Ax a 2 


A 2 a 2 


a/a 


Ca 


exp(— r/a) 


1 


ln2-| 




e -3/4 


8{a-r) 


2 


1 

2 


e -V2 


2e"^+ 1 /4 


(a/r) exp(— r/a) 


1 


In 2 


e"T 


1 


exp(-r 2 /a 2 ) 


2 


In 2 


e -7/2 


V2e-^ 2 



there are two different integrals, one of which being nonlocal. In three dimensions no such 
general closed form for the energy is available even near the threshold A = A c . 

The double integral in (|29|) resembles the one encountered earlier in the e-expansion for 
A c (e). For the spherically symmetric potential at small e one has 

A c (e) ~ A ie + A 2 e 2 + -- - , (32) 

where 

Al = ~J Q °°drrV(r) ' (33) 



1 / °° drr / °° dr'r' V(r) V(r') In ^ 

2 [J™drrV(r)} 3 



— " „ T ^: r< (34) 



( r > ( r <) is the greater (lesser) of r, r'). Then, after some straightforward calculations, we 
obtain another expression for C 

C = - exp [ In 2 - 7 - ^ ] , (35) 
a \ Ai/ 

where a is the mean range of the potential, defined by 

_ f™drr In rV(r) 

lna = — 7^1 — 777 \ • ( 36 ) 
j^drrV(r) 

We evaluate Ai, A 2 , a and C for several widely used potentials and the results are 
displayed in Table 1. Note that the values of the prefactor C are surprisingly simple. 

As mentioned earlier 0, equation (|32|) may be extrapolated to e = 1 to give a rather 
good estimate for the critical coupling in three dimensions. For the Yukawa potential this 
results e.g. in A c a 2 ~ 1 + ln 2 ~ 1.693 . . ., which is close to the exact result A c a 2 = 1.6798 . . .. 

It is also possible to use e-expansion in the same manner to evaluate the prefactor A 
in three dimensions. From ([3(]) it follows that A l ' e ~ C at small e. Extrapolating this 
result to e = 1 we see that actually C may be treated as a first approximation for A in 
three dimensions. For the square well potential this approximation gives Aa ~ 2e~ 7+1 / 4 ~ 
1.442 . . . the exact value being equal to Aa = tt 2 /8 = 1.234 . . .. While qualitatively correct, 
this first approximation is not very accurate. 
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Equation ( p29|) is valid for spherically nonsymmetric potentials as well, provided the range 
of the wavefunction k _1 is much larger, than the radius a of a potential. Hence equation 
( p9|) may be used in the problem of several attractive centers. 

In the case of two identical centers separated by a distance R we may write 

V(t) =v(r)+v(r + R) , (37) 



and upon substituting this potential in fl3l|) we obtain that C ~ exp(— = lni?) at R 3> a for 



arbitrary short-range v(r). Then for the energy Ei of a weakly bound state at a <C R <C k 1 
equations (p9|) and (|3T| ) yield 

£ 2 = -2e"^, (38) 

where kq is the square root of the binding energy on one center (given by (|29| ) and ( |3~1~D with 
the replacement V(r) — ►■u(r)). This is obviously the energy of the symmetric state. Note, 
that the energy E2 depends on the details of the interaction only through kq. 

It is interesting that the effective long range force between two centers, resulting from 
the bound state is of the Coulomb type. In three dimensions the corresponding energy is 
known to behave as 1/R 2 This is related to the collapse of three particle system with 
zero-range interaction, known as Thomas effect fL9f . Less singular behaviour of the energy 
(P8|) at small R is related to the absence of the Thomas effect in two dimensions [M . For 
separable potential of a particular type the dependence Kq/R in two dimensions was derived 
in JT2]. We see now that this formula is quite universal. 
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APPENDIX: A 



For the spherically symmetric potential the radial part of the 2 + e dimensional 
Schrodinger equation reads 

{-^t/ +1 Tt + xv{r) ) * (r) = mr) (A1) 

The substitution \1/ = <fir~( 1+e ^ 2 puts equation ( |A1|) in the form 

{^-^--mr)+Ey^0 (A2) 

This is obviously the radial equation for the wave function with non-zero angular momentum 
I in three dimensions with 1(1 + 1) = (e 2 — l)/4, i.e. I = (e — l)/2. 
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Next we proceed to the evaluation of the free particle Green's function, Ge{t), in 2 + e 
dimensions. This function certainly can be found elsewhere in the literature but, for the 
sake of completeness we give here its short derivation. The Green's function satisfies the 
equation 

(-A - E)G E {v) = 5{r) (A3) 
At r ^ we have for / = r^G^r) 

( — + 1 - - f2 +e] f - (A4) 
\ dr 2 r dr 4r 2 J 

Then, at E < 0, /(r) (which goes to zero as r — > oo) is proportional to the modified Bessel 
function K e /2(nr), where K = y/—E, hence 

G E (r) ~ r-^ 2 K e/2 ( K r) . (A5) 

At E — the Green's function reduces to the fundamental solution of the Laplace equation, 
which in 2 + e dimensions is given by equation (§). At nr — > and e < 2 we have 

K e/2 (.r) = 1 ^^j-y^ = (A6) 
1 K ' 2 sin(e7r/2) v ; 
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2sin(e7r/2) 



(^/ 2 )" £/2 _ (^/ 2 ) £/2 , 0(Ur?-^) 

r(i-e/2) r(i + e /2) +c;((Krj j 



Using this asymptotics, identity r(z)r(l — z) = it/ sm(irz) and comparing (|A5[) with (§|) we 
can fix unknown constant in ( |A~5|) 

1 

(27T)l+^2 V 



GE(r) = 7WZ ^(-Y /2 K e/2 (nr). (A7) 



Small k behaviour of ([AT]) gives rise to the expansion ( |l~4"D for (5Gb = G# — Go in the main 
text. 



APPENDIX: B 



To derive the separable potential flIB|) one can start from the following separable decom- 
position of a local potential V(r) 

V(t) = j £ <r n V\i> n ){,l> n \V , (Bl) 

n 

where the set of functions \ip n ) is determined from the eigenvalue equation 

(A + E)~ l Vip n = r] n (E)ip n , (B2) 
and a n = ±1 depending on the sign of (ip n \V\i[} n ) (see e.g. PTJ). For negative E this equation 



is in fact the Schrodinger equation for the bound states, where ip n is the wave function of a 
bound state with energy E in the potential l/rj n (E) V(r), i.e. 
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(-A + l/ Vn V(r))ij n = Eij n . (B3) 

From this equation one can derive a normalizing condition, which in this case is known to 
be (ip n \V(r)\ip m ) = a n 8 nm 0. 

Since we are interested in the weakly bound state with E ~ 0, we may now take the limit 
E — > in equation (|B3|), i.e. we may define the set ip n with respect to the zero energy. In 



this case ^ n 's are the threshold wavefunctions and the corresponding values of l/7/ n (0) are 
critical values of coupling constant. The wavefunction ipo of the first ground state at E — > 
obeys the zero energy equation (||), where A c = l/r)o (0) and rjo is the largest eigenvalue in 
QB2|). We do not know ?/v's exactly for an arbitrary potential, but e.g. for the attractive 



square well of radius a in three dimensions one can easily obtain for the s-states 



- sin \f\^r r < a 
r (_l)ni r>a 



(B4) 

A n = l/r} n = (±7r + irn) 2 ± . 

In the critical region A ~ A c the wavefunction of the first bound state is very close to 
tpo, so it seems natural to retain only the term with ipo in the expansion (|B~T| ) as a zero 
approximation. This approach is similar to the pole approximation in the scattering theory, 
valid for the resonance scattering when there exist a weakly bound state. 

Assume next that (V'ol^l'^o) < so that cr = —1. Then, in the vicinity of the critical 
point, where the transition from zero to one bound state occurs, we arrive at the separable 
potential (111 



V ~ V sep = - V\i> ) (*Po\V , (B5) 
where i/jq is normalized by the condition ([|) . 
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